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ABSTRACT 


An  exact  expression  for  the  total  intensity  of  the  pure 
rotational  band  system  of  a  rigid  symmetric  top  molecule  is 
derived.  An  approximate  expression  for  the  mean  spectral 
absorption  coefficient  is  obtained  and  is  compared  with 
previously  published  results.  The  first-order  effects  of 
non-rigidity  and  mechanical  anhaimonicity  on  the  band  intensity 
are  investigated  in  the  special  case  of  a  diatomic  molecule 
and  found  to  be  small  in  general. 
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INTRODUCTION 


The  total  intensity  of  the  pure  rotational  band  system  of  a  molecule 
is  defined  as  the  sum  of  the  intensities  of  all  transitions  of  the  molecule 
in  which  the  vibrational  quantum  number  does  not  change.  We  will  start 
with  the  expression  for  the  intensity  of  a  single  rotational  transition 
of  a  symmetric  top  molecule  and  sum  this  expression  over  all  the  quantum 
numbers  involved  to  obtain  an  exact  expression  for  the  total  band  intensity. 
An  approximate  expression  will  be  derived  for  the  mean  absorption  coefficient 
for  which  the  exact  summations  involved  cannot  be  performed. 

For  clarity  of  presentation,  the  following  development  neglects  the 
effect  <-f  alternating  statistical  weights,  and  strictly  applies  only  to 
an  accidentally  symmetric  top  molecule.  However,  the  proof  for  the 
general  case  of  a  symmetric  top  molecule  of  point  group  C^v  (such  as  the 
C^v  molecules  -  NH^,  CH^Ci,  etc.)  proceeds  in  exactly  the  same  manner, 
but  requires  more  cumbersome  equations  and  additional  definitions. 

TOTAL  BAND  INTENSITY 

The  intensity  of  a  pure  rotational  line  of  a  symmetric  or  asymmetric 
top  is  given  by^ 


qvjVm' 

vJtm' 


8tt^  ^vJtM  2 
3hc  p  ^o 


»[  1-exp (-  ||  w)]  £  |ft 


Fg1  JtM,J  t  M 


(1) 


where 

|  exp[  -E  (v,  J,  T  ,M)/kT]/ Q 


(2) 


is  the  number  of  molecules  per  unit  volume  and  unit  pressure  in  the 
state  described  by  the  usual  quantum  numbers  (v,J,t,M), 
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Q=  £  exp[-E(v;J,T,M)/kT] 


vJtM 


is  the  total  partition  function,  and  the  are  the  direction  cosines 
between  the  space -fixed  F  and  rotating  g  axes. 

If  Eq.  (l)  is  summed  over  M  and  M/,  we  obtain 


vjV  _  8rt^  N  2  e> 
vJt  “  3hc  p  ^o 


iliMl  «[l-e*p(-  »)]  ^  l»Pg|2jTM,j'TV,  (4) 


FM4' 


since  the  energy  levels  and  frequencies  are  independent  of  M  and  M  . 
If  we  consider  the  top  rigid,  i.e., 

E(v,J,t)  =  hc[G(v)  +  F(J,t )]  , 


and  sum  over  v,  we  obtain 


_/  / 
T 


-  I  “o  Cxp[  ^  -(-J  —  U)[l-exp(-  H  •)]  £  |*Fg|2  ,  (6) 


where  is  the  rotational  partition  function 


%  -  I  e*p[-  If  • 


We  now  specialize  to  the  case  of  the  symmetric  top,  for  which  we 


have 


F(J,K)  =  BJ(J+l)  +  (A-B)  K2, 


where  the  symmetric  top  quantum  number  K  (which  runs  from  0  to  J)  is 
now  used  in  place  of  t.  For  a  J  -*  J+l  transition,  we  have  (since  AK  =  0) 
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u>  =  2B(j+l), 


(9) 


y  i«  i2  ,  -  (p.6  j  lj+.d2-*5-  ,  (10) 

L  1  Fg'  JKM,  J+1K  M  "  °0K;  J+l 

fmm' 

From  Eq.  (6),  we  obtain 

„J+1  K  8n3  K  .2  eXP(-  itBJ(J+l)  +  <*’B> 

J  K  =  3hc  p  ^o 

X  2B(J+1)  {l-exp[ -  g§  2B(J+l)]|  (2-6OK).  (11) 

Hie  quantity  Sj+1  is  now  obtained  by  summation  over  K: 

sj+1  ■  is  f  4 1  {exp[-  t§  BJ<J+1)]  -exp[-  §  b(j+i)(m)]} 

J 

X  1  (2-5ok)  [(J+lf-K2]  exp[-  |§  (A*B)  K2] .  (12) 

K=0 

Hie  summations  over  K  in  Eq.  (12)  cannot  be  expressed  in  closed 
form.  We  will  approximate  the  sums  later  by  integrals,  and  derive  an 
approximate  expression  for  the  absorption  coefficient.  However,  we 
will  continue  the  present  development  in  order  to  obtain  an  (exact) 
expression  for  tne  band  intensity. 

If  we  define  the  quantities  f(j)  and  g(j)  as  follows: 

J 

f(j)  =  £  (2-6qk)  exp[-  |§  (A-B)  K2]  (13) 

K=0 


h 


and 


g(j)  =2^1^  exp[-  ||  (A-B)  K2], 


(1U) 


K=0 


then  we  have 


3 

sj+1  “  3?  p  iexp["  lfBJ(J+1>j  -exp[-  ^(j+i)(J+2)]}[(J+i)2f(J)‘6(J).. 


and  Eq.  (7)  becomes 

00 

Or  =  £  (2J+l)  exp[-  |^J(J+1)]  f(j).  (16) 

J=0 

We  note  here  that  if  the  derivation  is  carried  through  for  a  non- 

accidentally  symmetric  top  molecule,  an  additional  nuclear  spin-dependent 

factor  in  the  statistical  weights  must  be  considered.  (This  factor  has 

implicitly  been  assumed  constant  here  and,  hence,  has  been  neglected.) 

2 

This  factor  is  dependent  only  on  the  quantum  number  K.  Extended 
definitions  analogous  to  Eqs.  (13 )  and  (l4)  can  then  be  made  and  carried 
through  in  the  same  manner  to  produce  the  same  results. 

The  total  rotational  band  intensity,  given  by 

GO 

«(T)  =  l  Sj+1,  (17) 

J=0 

is  obtained  by  substituting  Eq.  (15)  into  Eq.  (17)  and  rearranging  terms: 

■w  ■  1 1  ») 
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-L  +  Iexp[-hcBj(J^l)/kT]{j2[f(j)-f(J-l)1-g(j)-t-g(j-l)} 


I 


exp[-hcBJ(J+l)/kT]  (2J+l)  f(j) 


(15) 


.  (18) 


But  from  the  definitions  of  f(j)  and  g(j),  we  have 


=  2J2  exp[-  g  (A-B)  J2]  (19) 


and 


[f(j)-g(j-l)]  =  2J2  exp[-  jg  (A-B)  J2] 


(20) 


Thus  the  curly-bracketed  term  in  Eq.  (18)  vanishes  identically,  leaving 


«<T)  =  #  2  /  B. 
'  3hc  p  o 


(21) 


Since  — 
P 


kT 


,  this  can  be  written 


<y(T)  = 


O  O 

lfirr-3  p,  B 
o 

3hc~kT 


(22) 


or 

a(273°K)  =22.?  p,2  B  cm" 2atm-1,  (23) 

where  B  is  in  cm  1  and  p,Q  is  in  debye  units.  Thus  we  find  that 

\  a(Tx)  =  T?  a(T2),  (24) 

that  is,  the  (density-corrected)  intensity  of  a  pure  rotational  bane 
system  of  a  rigid  symmetric  top  molecule  is  independent  of  temperature . 

We  note  here  that  the  expression  for  the  total  band  intensity  is 
a  function  of  B  but  not  of  A.  This  is  related,  at  least  in  part,  to  the 
fact  that  the  frequency  of  the  radiation  depends  only  on  B.  The  limiting 
case  A-*®  was  treated  by  Golden,  who  determined  the  absorption  coefficient 
as  a  function  of  frequency  for  the  rigid  rotator  model.  Integration  of 
Golden’s  Eq.  (5)J  over  frequency  yields  an  expression  for  total  band 
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intensity  which  is  consistent  with  Eq.  (22)  within  the  approximation 
involved  in  replacing  the  summation  over  J  by  an  integration. 

APPROXIMATE  MEAN  ABSORPTION  COEFFICIENT 
The  mean  absorption  coefficient  k  is  given  by  Sj+  /d,  where  d  -  2B 
is  the  line  spacing,  so  that  we  have  from  Eq.  (15 )  and  Eq.  (9) 

k  =  s^/a  =  ^  {«*p[-  $p(J.l)]  -exp[-  |^(J+l)(J+2)]} 


X  [(J+l)2f(J)-g(J)]. 


The  sums  involved  in  the  definitions  of  f(j),  g(J)  and  cannot  be 
performed  exactly.  If  we  approximate  the  sums  by  integrals,  we  obtain 

M) 

f(j)  *  J  exp[-  ||  (A-B)  K2]  dK 

-(J+i) 


^  ■  -■  erf L /  ||  (A-B)  (J+&)]  , 

m  a-b) 


and,  similarly, 


2l||(A-B)]' 


ev/J  |§(A-B)  (J+i)]‘  exp['  |(a-B)(J^)2J-  (2V) 


i(A-B) 


To  the  same  degree  of  approximation,  we  have  for  the  rotational 


partition  function 


(j+£ ) 


%  -  |  (2J+1)  exp[-  ||bJ(J+1)]  |  «p[-  jgtA-BjK2]  BKdJ 

4  -(J+^) 


%  ■  =x4  §•)  7==  I  2Iexf[- 1 BI®] 

m*- b) 


evty  g(A-B)lJ  dl. 


r®  2  2  * 

But  since  x  exp(-x  )  erf (ax)dx  =  ^  a (a  +l)  ,  we  have 

v  _ 


kTv3^2  1 


if  we  set  exp(^  ^  B)  «  1.  Thus  the  use  of  this  customary  expression 
for  Qp  is  completely  consistent  with  the  assumed  model  and  approximations; 
further  comment  will  he  made  on  this  point  in  discussing  the  limiting 
case  A-^. 

Thus  we  obtain  the  following  explicit  expression  for  k: 


k  -  (|?  5  4)  ^#3/2  |^(^)]-exp[-  |§B(J+l)(J+2)]} 


tf(A-B) 


(J+l)2  erf (/ ||(A-B) '  (J+4)] 


2[^(A-B)]' 


erf[./^(A-B)  (J+4)]  + 


exp[-  ||(A-B)(j+4)2j},  (30) 


or,  to  a  degree  of  accuracy  consistent  with  the  previous  approximations. 


k  ■  (|£  |  4)  /SW-  ||BJ(j+l)]-exp[-  j§B(J+l)(J+2)]} 


x  (J+l)2  G[/(A-B)  £§  (J+l)]  , 


where 
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G(z)  =  Jn  (l-  erf(z)  +  -  exp(-z^). 
2z  z 


(32) 


The  power  series  representation  of  G  is  given  by 

00 

- 4  I  a(S5fi)C*3T %ZaA’  (33) 

n=0 

Equation  (31)  is  not  consistent  with  the  expression  for  the  absorption 
coefficient  presented  by  Fenner^  [Eq.  (11-134)^].  There  appear  to  be  two 
causes  of  discrepancy. 

First,  the  induced  emission  term  [l-exp(-hcu>/kT)]  was  implicitly 
assumed  to  be  approximately  unity  in  the  derivation  of  Eq.  (11-134)^. 

Such  an  approximation  is  not  valid  in  gene,  xl  for  a  pure  rotational  band 
as  can  be  shown  by  the  following  argument:  Consider  a  spherical  top 
molecule  (for  simplicity).  The  population  of  a  given  J  state  is  approxi¬ 
mately  proportional  to  J2  exp(-gJ2),  which  is  a  maximum  for  3J2  =  1.  For 
this  value  of  J,  the  induced  emission  term  [l-exp(-23  J)]  is  equal  to 
[l-exp(-2/j)] .  Even  at  temperatures  low  enough  for  the  maximum  population 
to  lie  in  the  J  =  3  *tate,  the  induced  emission  term  is  equal  to 
[l-exp(-2/3)]  rs  0.49.  With  increasing  temperature,  the  induced  emission 
term  decreases;  thus  approximating  this  term  by  unity  is,  in  general,  not 
reasonable. 

Secondly,  Eq.  (31)  differs  from  Bi.  (11-134)^  by  a  factor  of 
approximately  2  or  3*  This  appears  to  result  from  the  use  of  a  summation 
rule  [Eq.  (7.116a)  ]  in  which  a  sum  of  matrix  elements  is  taken  over  the 
upper  state  values  of  J-l,  J,  and  J+l,  by  analogy  to  the  case  of  a 
vibration- rotation  band  system.  In  such  a  system,  a  lower  state  J  can  go 
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to  an  upper  state  J-l,  J,  or  J+l.  However,  in  the  pure  rotational  band 
system  under  consideration,  the  only  corresponding  transition  is  J  -*  J+l. 

o 

In  the  limit  of  A-+“,  Eq.  (3l)  is  consistent  with  Golden's  Eq.  (5), 
to  the  degree  of  approximation  previously  noted.  Since  it  can  be  seen 
that  the  approximate  expressions  for  f(J),  g(j)  and  become  arbitrarily 
poor  for  large  A,  we  will  comment  briefly  on  this  point. 

We  note  that  in  the  limit  as  A-*®,  we  have  from  Eqs .  (13) ,  (1*0,  an<^ 

(16) 

f(j)  -  1  ,  (3*0 

g(J)  -  0  ,  (35) 

and 

CO 

Or  -  £  (2 J+l)  exp[-  — BJ(J+l)],  (36) 

J=0 


i.e.,  the  rotational  partition  function  for  a  linear  molecule.  Approximating 
this  sum  by  an  integral  in  the  usual  way  yields 


(37) 


Thus,  in  this  limit  we  have 

f(j)/QR  -  hcB/kT 


(38) 


and 

g(j)/Qp  -  0.  (39) 

But  exactly  the  same  results  are  obtained  by  taking  the  ratios  of  the 
approximate  expressions  used  (i.e.,  the  ratios  of  Eqs.  (26)  and  (27)  to 
Eq.  (29)).  Only  these  ratios  are  used  in  the  present  work;  thus  the 
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limiting  case  of  A r*°  is  a  valid  limit  although  the  separate  expressions 
for  f(J),  g(j)  and  becane  arbitrarily  poor  approximations  for  large  A. 

EFFECTS  OF  CENTRIFUGAL  STRETCHING 

Wf  will  consider  here  only  the  special  case  of  a  diatomic  molecule; 
a  general  formulation  for  polyatomic  molecules  is  not  feasible  since 
each  different  molecular  configuration  presents  a  different  geometrical 
problem.  It  is  felt  that  the  magnitudes  of  the  effects  determined  for 
the  diatomic  molecule  will  be  representative  at  least  of  those  of  the 
simpler  polyatomic  molecules. 

We  will  consider  the  effect  of  nonrigidity  of  the  molecule  on  the 
band  intensity  through  (a)  its  effect  on  the  energy  level  structure 
(and,  hence,  on  the  emitted  frequencies)  and  (separately)  through  (b) 
its  effect  on  the  electric  dipole  moment. 

One  would  expect  that  since  centrifugal  stretching  tends  to  lower 
the  energy  levels  and  decrease  the  net  effect  of  (a)  would  be  to 

lower  the  total  band  intensity.  Also,  since  the  stretching  may  tend  to 
increase  the  dipole  moment  (at  least,  for  diatomic  molecules,  in  many 
cases),  one  would  also  expect  that  the  net  effect  of  (b)  would  be  to 
increase  the  total  band  intensity.  These  two  effects  will  be  investigated 
quantitatively . 

The  energy  levels  of  a  diatomic  molecule  (with  the  first  approximation 
to  the  effect  of  centrifugal  stretching)  are  given  by^ 

F(J)  =  BJ(J+1)-DJ2(J+1)2.  (40) 

The  second  term  is  assumed  to  be  a  small  correction.  For  this  simple 
model,  D  is  related  to  the  rotational  and  vibrational  constants:^ 
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(41) 


D  /2B\  2 

B  =  (~)  =  Y  * 

v 


The  frequency  of  radiation  is  then 


u,JT+1  =  F(J+1)-F(J)  =  2B(J+l)-4D(J+l)3. 


(42) 


We  then  have  for  the  intensity  of  a  J  —  J+l  transition 

3  2 

SJ+1  *  [ihT  S  ^  C2B(J+1)-4D(Jt1)3]  <J+1> 

x  jexp[ -0 J ( J+l)]exp[ 6  ( J+l)^l -exp[ -0 ( J+l) ( J+2 )]exp[ 6 ( J+l)^ ( J+2 )2]  j  (43 ) 

(using  the  notation  0  =  jj|^B  and  6  =  j^D). 

If  we  expand  the  exponential  functions  involving  D,  retaining  the 
first-order  terms,  and  proceed  as  before,  we  obtain  for  the  band  intensity 

r  cj+i 

“  =  LSJ  : 

-QnV 

o 


-  8D^  (2 J+l) J ( J+l)exp[ -0 J ( J+l)] 

+  2B6^  (2J+!)J2(J+I)2exp[-0J(J+I)]|  .  (44) 


To  the  same  approximation  we  have 


Qj,  =  V  (2J.-l)exp[-SJ(J+l)][l+6J2(J+l)2] 


=  "S' 


where 


=  0)  =  0"1.  Thus  we  have 


(45) 
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I,lii1||l|||||i|  i-mumifiltii-  ^ WWiii ii 'H< . .  '“’"'"'I . utf . .  . . . . . . . . . . . . . . 


-2- 


%  =  Q^[ 1+263  ""]. 

In  terms  of  Q^,  we  have  for  the  band  intensity 


0^(1+263  ) 


(46) 


(47) 


from  which  we  obtain 


a  =  V%^  ~ 1  - ^—p  [2B-8D3-142B63"2] 

L  3hc  pJ  li0,a-2  L  J 


14263 


(46) 


or 


-8n3u2 


(49) 


We  now  consider  the  effect  of  the  varying  dipole  moment.  We  write 


d  =  4  Ar, 


(50) 


and  express  Ar  in  terms  of  the  molecular  constants : 


hr  =  ^  r  J(J41)  =  y2reJ(j4l). 

% 


(51) 


Thus  we  find 


rp,  p 

M-  =  dQ[l+  “ 1  V  J(J4l)] 


(52) 


or 


2  2 


d  = 


d^  [1+2  J(J+l)]  , 


(53) 


where 


e  =  ivSv 


(54) 
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o  2 

If  this  value  of  k  is  substituted  for  uo  in  Eq.  (ll),  we  obtain 
after  similar  manipulation, 


*  =[ 


q_  3  2 

™  ^o  m 


3hc  pJ 


2Bi  1  +  4 


•£_  _kin 
e  hcBJ  * 


(55) 


If  we  combine  Eqs .  (49)  and  (55),  we  find,  for  the  net  effect  of 
centrifugal  stretching, 

3  2 

.  S]  S  k  D  JcT  4  •£_  _kTl 

L  3hc  pJ  ■!_  B  hcB  8  hcBJ 


(56) 


D  2B  2 

If  we  assume  further  that  —  =  ( — )  (=y  )  [Eq.  (4l)],  we  have 

L  o>v 


.3.2 


r^o  tnoJ\  ,,  2  , .  lv  kTl 


(57) 


2  ^[p  -5 

Investigating  orders  of  magnitude,  we  find  that  since  4y  ~  T  x  10  ^ 

3 

for  most  diatomic  molecules,  the  correction  is  small  for  T~  10  .  We 
note  that  (l-  ^)  is  not  large  unless  8  »0,  which  implies  that  uQ  (and  cy) 
are  small.  Although  then  1(1-  hi  *  I!  1,  there  is  a  multiplicative 

D  O 

2 

factor  of  uo  in  Eq.  (57),  which  will  still  dominate.  Hence,  for  9  ^  0, 
the  percent  correction  my  be  large,  but  it  is  small  on  an  absolute  basis. 

For  the  other  special  case  of  9  =  +1  (which  is  approximately  true  for 
HC£  and  HF),  the  net  effec  of  centrifugal  stretching  vanishes. 

EFFECTS  OF  MECHANICAL  ANHARMONICITY 
We  next  consider  the  fact  that  the  rotational  constants  my  vary 
with  the  vibrational  quantum  number  (as  a  result  of  the  mechanical 
anharmonicity  of  a  real  molecule). 

If,  in  the  previous  development,  the  rotational  constants  are  not 
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considered  independent  of  v  as  in  Eqs .  (5)  and  (8),  and,  instead,  the 
expression 


B(v)  =  B  -  a  v 
x  7  o  e 


(58) 


is  used,  we  obtain  (after  a  similar  sequence  of  manipulations  and 
approximations ) 


a 


L  3hc 


(59) 


In  effect,  the  rotational  constant  B  is  replaced  by  the  (population- 
weighted)  average  value  of  B(v)  at  temperature  T.  The  correction  term 
is  of  the  order  of  1-.0002  for  T  ~  lCr  . 

In  order  to  obtain  the  corresponding  correction  for  the  dipole 
moment,  we  note  that  in  a  vibrational  state  defined  by  the  quantum 
number  v,  we  have 

^av  ■  <[Mb  +  “l  ir]a>av 


o  M>  -i 

\i  (1+2  —  <Ar)  ). 
*0  '  u  av7 


(6o) 


-2v 


To  relate  (Ar)  to  v,  we  note  that  since  B  «  (r  )  ,  we  have 

av  h.  v 


av" 


B(v)  =  B  -  a  v  »  <r“!)  r  '2  (1  -  2(Ar>av/r  ). 


(61) 


Thus,  we  have,  approximately 

2<Ar)av/re  =  oev/BQ  (62) 

and 


IS 


(63) 


<4v  =  l 

o 


Proceeding  as  before,  we  find  that  we  have  (considering  only  the 
effect  of  the  dipole  moment  varying  with  each  vibrational  state ) 

r8TT^o  1TL_  L  lffe[  /he  x  .T1! 

a  -  Liet  pJ2Bo  l1  +  0  r  LeXp(kT  tUo)"1J  J  * 


(64) 


Oil  combining  the  effects  of  variable  B  and  p,,  we  obtain 


(65) 


The  comments  made  after  Eq.  (57)  also  apply  here,  although  the  magnitude 
of  the  correction  is  considerably  smaller,  in  general. 

Thus  the  expression  for  the  rotational  band  intensity  combining  the 
corrections  for  the  variation  of  B  and  p,  with  both  v  and  J  can  be  written 
as  follows : 


r8rT^p,2  .  / 


3hc 


2  kT 


hcB 


+  r  [exp(Ha,o)-1]’1)}-  (66) 
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